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Abstract 

We discuss the grazing collision limit of certain kinetic models of 
Bose-Einstein particles obtained from a suitable modification of the one- 
dimensional Kac caricature of a Maxwellian gas without cut-off. We re- 
cover in the limit a nonlinear Fokker-Planck equation which presents many 
similarities with the one introduced by Kaniadakis and Quarati in [T3|. 
In order to do so, we perform a study of the moments of the solution. 
Moreover, as is typical in Maxwell models, we make an essential use of 
the Fourier version of the equation. 
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1 Introduction 

The quantum dynamics of many body systems is often modelled by a nonlinear 
Boltzmann equation which exhibits a gas-particle-like collision behavior. The 
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application of quantum assumptions to molecular encounters leads to some di- 
vergences from the classical kinetic theory [5] and despite their formal analogies 
the Boltzmann equation for classical and quantum kinetic theory present very 
different features. The interest in the quantum framework of the Boltzmann 
equation has increased noticeably in the recent years. Although the quantum 
Boltzmann equation, for a single specie of particles, is valid for a gas of fermions 
as well as for a gas of bosons, blow up of the solution in finite time may occur 
only in the latter case. As a consequence the quantum Boltzmann equation 
for a gas of bosons represents the most challenging case both mathematically 
and numerically. In particular this equation has been successfully used for com- 
puting non-equilibrium situations where Bose-Einstein condensate occurs. From 
Chapman and Cowling [5j one can learn that the Boltzmann Bose-Einstein equa- 
tion (BBE) is established by imposing that, when the mean distance between 
neighboring molecules is comparable with the size of the quantum wave fields 
with which molecules are surrounded, a state of congestion results. For a gas 
composed of Bose-Einstein identical particles, according to quantum theory, the 
presence of a like particle in the velocity-range dv increases the probability that 
a particle will enter that range; the presence of f{v)dv particles per unit volume 
increases this probability in the ratio 1 -I- 6f{v). This fundamental assumption 
yields the Boltzmann Bose-Einstein equation 

df 

-^ = QQBE{f){t,v), tER+, veR\ (1) 

where 

QQBE{f){t,v)^ f B{v-v,,u){ff:{l + Sf){l + 5f,) 

JR3XS2 (2) 

-fMl + Sf){l + Sf:))dv,doj, 

where as usual we denoted 

f^fiv), = .f' = fiv'), fi = fK), 

and the pairs {v,v,,) (respectively {v',vl)) are the post- (respectively pre-) col- 
lision velocities in a elastic binary collision. In Q B{z, u) is the collision kernel 
which is a nonnegative Borel function of |z|, | < z,uj > \ only 

B{z,u:)^B(^\z\,^^^y {z,uj)eR''xS'. (3) 

The solutions f{v,t) are velocity distribution functions (i.e., the density func- 
tions of particle number), 5 = {h/m)^/g, h is the Planck's constant, m and g 
are the mass and the "statistical weight" of a particle (see [17J for details). 

For a non relativistic particle, by setting v{p) = p/m, the collision operator 
Qqbe can be rewritten in general form as follows |22l 123] 

QQBE{f){t,p)= f w{p,p,,p',p:){ff:ii+sf){i+sf,) 

JR9 (4) 

- fM^ + Sf'){i + Sfl))dp,dp'dp: 

where is a nonnegative measure called transition rate, which is of the form 

wip,p,,p\p:) = nip,p,,p\p:)sip,+p.,-p' -p:)di£{p)+£(p.,)-£ip')-£ip:)), 
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where 5 represents the Dirac measure and £{p) is the energy of the particle. The 
quantity Wdp'dp'^ is the probabihty for the initial state (j),p*) to scatter and 
become a final state of two particles whose momenta lie in a small region dp'dp'^ . 
The function 51 is directly related to the differential cross section (see ([3])), a 
quantity that is intrinsic to the colliding particles and the kind of interaction 
between them. The collision operator ^ is simplified by assuming a boson 
distribution which only depends on the total energy e = £{p). In this last case 
/ = /(e, t) is the boson density in energy space. 

Together with the Boltzmann description given by the collision operators 
(El)-®; other kinetic models for Bose- Einstein particles have been introduced 
so far. In particular, a related model described by means of Fokker-Planck type 
non linear operators has been proposed by Kompaneets [TS] to describe the 
evolution of the radiation distribution f{x,t) in a homogeneous plasma when 
radiation interacts with matter via Conipton scattering 



In that context the coordinate p represents a momentum coordinate, p = \p\. 
More precisely, an equation which includes ([S]) as a particular case is obtained 
in |15| as a leading term for the corresponding Boltzmann equation under the 
crucial assumption that the scattering cross section is of the classical Thomson 
type (see [TU] for details). 

The fundamental assumption which leads to the correction in the Boltzmann 
collision operator namely the fact that the presence of f{v)dv particles per 
unit volume increases the probability that a particle will enter the velocity range 
dv in the ratio 1 + Sf{v), has been recently used by Kaniadakis and Quarati 
[T51 [T^ to propose a correction to the drift term of the Fokker-Planck equation 
in presence of quantum indistinguishable particles, bosons or fermions. In their 
model, the collision operator ^ is substituted by 



Maybe the most remarkable difference between the kinetic operators ^ and 
([S]) is that, while the former is such that mass, momentum and energy are colli- 
sion invariant, the latter does not admit the energy as collision invariant. This 
suggests that the operator ^ would not result directly through an asymptotic 
procedure from the Bose-Einstein collision operator ^ , but instead from some 
linearized version, where only the mass is preserved under the collision mecha- 
nism. 

For a mathematical analysis of the quantum Boltzmann equation in the 
space homogeneous isotropic case we refer to |S1 El UHl HZl UH] ■ We remark that 
already the issue of giving mathematical sense to the collision operator is highly 
nontrivial (particularly if positive measure solutions are allowed, as required 
by a careful analysis of the equilibrium states). All the mathematical results, 
however, require very strong cut-off assumptions on the cross-section |10[ I17| . 

Also, accurate numerical discretization of the quantum Boltzmann equation, 
which maintain the basic analytical and physical features of the continuous 
problem, namely, mass and energy conservation, entropy growth and equilibrium 
distributions have been introduced recently in [21 HF] . Related works [111 HOI [S] 
in which fast methods for Boltzmann equations were derived using different 





(5) 



Qj.p(/)(t,«) - V - [V/ + ^;/(l + (5/)]. 



(6) 
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techniques like multipole methods, multigrid methods and spectral methods, 
are relevant to quote. 

At the Fokker-Planck level, the qualitative analysis of the Kompaneets equa- 
tion described by the operator ^ has been exhaustively studied in [2, while 
the numerical simulation has been done by Buet and Cordier [3] • To our knowl- 
edge, the mathematical study of the Fokker-Planck equation ^ introduced by 
Kaniadakis and Quarati |13| has been done only very recently [4J, where the 
one-dimensional version of ^ has been studied. 

In the case of the quantum Boltzmann equation the asymptotic equivalence 
between the binary collision operators Q and the Fokker-Planck type op- 
erators ([5]) and ([6]) is unknown. This is not the case for the classical binary 
collisions in a elastic gas, where the asymptotic equivalence between the Boltz- 
mann and the Fokker-Planck-Landau equations has been proven rigorously in 
a series of papers by Villani |25l r26j by means of the so-called grazing collision 
asymptotics. 

The same asymptotic procedure, in the case of the one-dimensional Kac 
equation |TT] , showed the asymptotic equivalence between Kac collision operator 
and the linear Fokker-Planck operator The method of proof in [24] takes 
advantage from the relatively simple structure of Kac equation. Taking this into 
account, in order to establish the asymptotic connection between the Boltzmann 
equation for Bose Einstein particles and its Fokker-Planck description, we will 
introduce a one-dimensional kinetic model in the spirit of Kac caricature of a 
Maxwell gas with a singular kernel. Then we will study the grazing collision 
limit of the equation, which leads to a Fokker-Planck type equation in which the 
drift is of the form of equation ([5]), but the coefRcient of the (linear) diffusion 
term depends on time through the density function. More precisely, the Fokker- 
Planck collision operator reads 

Qppim, v) = + + Sf)), (7) 

where 

Mf)= I v^f{v,t){l + Sf(v,t))dv. (8) 
Jm 

The paper is organized as follows. In the next Section we will introduce the 
model, together with some simplifications (mollified model) that allow to prove 
in Section |3] existence of a weak solution. Then in Section 21 we will focus on 
the moments of the solution and on some regularity properties. In Section [SI 
we will deal with the grazing collision limit. The last part of the work. Section 
[6] will contain some results for the non mollified model. 



2 The Kac caricature of a Bose Einstein gas 

The simplest one-dimensional model which maintains almost all physical prop- 
erties of the Boltzmann equation for a Bose-Einstein gas can be obtained by 
generalizing Kac caricature of a Maxwell gas to Bose-Einstein particles. This 
one-dimensional model reads as follows 

df 

-^ = QQBE{f){t,v), teR+, veR, 
f{0,v) = fo{v), 
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where 



QQBEim,v) ^ J j{e) J^{f'fUi + Sf)ii + 6U) ^^^^ 

-ff.a+sn{i+sf:))dv,d9. 

For the sake of brevity, we used the usual notations 

f = f{t,v), f' = f{t,v'), u = f{t,v,), f: = f{t,v',), 

The initial datum /o is a nonnegative measurable function. The pre-collision 
velocities are defined by the Kac rotation rule [11], which is given by 

I w' = w cos 6' — sin ^^^^ 
cos 0. 



Collisions pip imply the conservation of the energy 



2,2 12 , 12 



(12) 



Let us observe that the system ((TT|) can be reversed and that we can write the 
post-collision velocities with respect to the pre-collision ones 

I ti = u' cos 6* + sin 

(13) 



The parameter 5 in (|T0|) is a positive constant. The choice 5 — would 
lead us back to standard Kac model, whereas 5 negative would lead us to the 
Boltzmann-Fermi-Dirac equation, whose features are very different from that of 
the Boltzmann-Bose-Einstein equation. 

The cross-section I3{9) is a function defined over {—^,^). In the original Kac 
equation [Tl], /3(0) is assumed constant, which implies that collisions spread out 
uniformly with respect to the angle 9. Following Desvillettes [1], we will here 
assume that the cross-section is suitable to concentrate collisions on the grazing 
ones (these collisions are those that are neglected when the cut-off assumption 
is made) . This corresponds to satisfy one or more of the following hypotheses 

HI: I3{9) is a nonnegative even function. 

H2: I3{9) satisfies a non-cutoff assumption on the form 

m - 1^ when 9 ^ 0, (14) 



with I <v <2. That is. 



whereas 



/3(6l)|sine'|de' = 



P{9)\sin9\''+''d9 < + 



for aU e > 0. 
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H3: is zero near and ^, namely there exists an eo > such that 

V^?e(-|,-|+eo]n[|-eo,|). 

In the case in which the classical Kac equation is concerned, the asymptotic 
equivalence between the non cut-ofF Kac equation and the linear Fokker-Planck 
equation as collisions become grazing has been proven in |24| . Hence, the pas- 
sage to grazing collisions in (fTU)l . would give us the correct Fokker-Planck type 
operator which leads the initial density towards the Bose-Einstein distribution. 

Due to the symmetries of the kernel (fTO)) and to the microscopic conservation 
of the energy it can be easily shown, at least at a formal level, that the 
mass and the global energy of the solution are conserved 

fit,v)dv = / foiv)dv, 

and 

v'^f{t,v)dv= / v'^fo{v)dv, 



for alH > 0. Moreover, if 

H{f)^ J^(^^{l + Sf) log(l + Sf)-f log dv 
denotes the Bose-Einstein entropy, the time derivative of H{f) is given by 

^^^^ " i /I ^^^^ L ^ '5/)(l+5/*), 

fMl + Sf'){l + Sfi))dvdv,d9 



ith 



(a - 6) log(a/&), a > 0, 6 > 0; 
r(a,6) = <( +00, a>0, b = or a = 0, 6 > 0; (15) 

0, a = b = 0. 

Then, since D{f) > the solution f{t) to equation ^ satisfies formally an 
if-theorem. H{f{t)) is monotonically increasing unless f{t) coincides with the 
Bose-Einstein density fsE, defined by the relationship 

Ibe{v) 



l + 5fBE{v) 



ae-"" , (16) 



where a and b are positive constant chosen to satisfy the mass and energy 
conservation for /be- 

It can be easily verified by direct inspection that the fourth order term in 
(fTII|) cancels out from the collision integral, so that it can be rewritten as 



QQBE{f){t,v) = I _m i {fai+6f+5u) ^^^^ 

- fMi + 6f' + 6fi))dv,de. 
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In trying to give a rigorous signification to equation several difficulties 
arise. In fact, our non-cutoff cross-section /3(0) does not allow us neither to use 
the same change of variable as in |17j, nor to use the same weak formulation as 
in [3] . A sufficient condition to give a sense to the collision kernel would be that 
/ G L°°(K+ X R). It would even be enough that such a condition hold for the 
quantum part, that is for the / involved in the terms of the form 1 + Sf. To 
satisfy that condition, we modify the quantum part by smoothing it. Let -0 be 
a moUifier, that is 

1. ?A e C,°°(R) 

2. V > 

3. J^ip{v)dv = 1. 
Then, let 

f(t,v)^ / f{t,v-w)ilj{w)dw ^ f{t,.)*y'4>. 



The function / is regular in the velocity variable, and relies uniformly in all the 
i''(R) spaces (for 1 < p < +oo) since the L^(K.) norm of f{t,.) is constant. 
Moreover, f{t, .) is as close as we want to f{t, .) in all these norms, provided V' 
is well chosen, so that our new model is nothing but an approximation of ([9]) : 

df 

-^ = QQBE{f), teR+,veR ^^^^ 



fiO,v)^foiv) 



with 



J-f Jr (19) 

This approximation still formally preserves mass and energy, while maintaining 
the same nonlinearity of the original collision operator. It has to be remarked, 
however, that both the validity of the iJ-theorem and the explicit form of the 
steady solution are lost. Other approximations can be introduced, which do not 
exhibit this problem. Among others, the operator 

■ (1 + sm + if.xi + «/')(! + sfDdv.m. 

preserves mass and energy, satisfies the iJ-theorem and possesses the right 
steady state. Unlikely, the nonlinearity of (PU]) is difficult to handle for our 
purposes. 

The remaining of our work will concentrate on the study of (jl8p . In Section 
131 we will study the existence of a solution to this problem. Then in Section 21 
we will focus on the moments of the solution and on some regularity. In Section 
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[51 we will finally study the grazing collision limit. Our work will end in Section 
[Hlwith some partial result for the non mollified model ©. 

Let us end this Section with a few notations. The functional spaces that will 
be used in the following, apart from the usual Lebesgue spaces, are the weighted 
Lebesgue spaces, defined, for p > 0, by the norm 



llffllLi(E)= / {l + \vrMv)\dv. 

Jr 

We will also need some Sobolev spaces, defined for < s < 1 by the norm 

\\9\\HHR)^\\9\\h + \g\h 

where 

'^'- = 71 — — 

Our convention for the Fourier transform is the following: 



and the inverse Fourier transform is given by 

f{v) - ^ / /(Oe™«de 
We will sometimes use the notations 



fo{v)dv 

JR 

and 

3 Existence theorems 

The goal of this Section is to prove the existence of a solution to the problem 
. To start with, we first consider the case of a cross-section with some cutoff. 
We have 

Theorem 1. Let /o G L^iM.) be a nonnegative function. Assume that the cross- 
section satisfies j3 G L^[—^; ^) and (HI). Then there exists a unique solution 
f G L°°(R+; ijl-I^)) the problem U8\}. which is nonnegative, and preserves 
mass and energy. 



Proof For f,g,he i^(M), let 



f'9'JTM^ + Sh + 5k)-fg/ -^^^ 



I/oIIli 7 ' Vll/oll 



+ Sh' + Sh' 



dv.de + f[ I g][ + 2S\\ij\\L^ I h] I mdO. 
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Let 

if =|l/o|Ui(l + 25|l^|U=o||/o|Ui) r /3(0)d0. 
Let us consider the following problem 

^+Kf = P{f,fJ), teR+,v€R ^21) 
/(0,«) = /o. 

The operator P : (-^2)'^ ^ -^2 trilinear, and satisfies the inequality 
\\Pif.9.h)hi < Cp\\f\\^\\g\\^\\hh. yf,g,he Ll{R), 

with 

Cp^J m^e + 

Assume for the moment that K — 1. Thanks to a theorem of [T4], there exists 
some T > such that there exists a solution / e L°°{0, T; Ljl^^)) to the problem 
(|2ip. This solution can be written as a Wild sum, which reads 

+00 

where 

n+j2+i3=fc— 1 

and 

/fc=0 = /o- 

The numbers 6^ are the coefficients of the Taylor expansion of 

^ +00 

One can easily see that all the bk are positive. Moreover 

0<f,g,heLl =^ P{f,g,h)>0. 

Thus, the solution / of ([?!]) is nonnegative. Moreover, owing to the definition 
of P(-, •, •), one can verify that this solution preserves the mass. From that we 
deduce that / is solution of (fT8)) on (0,T). But it preserves the mass, and it 
relies in L\, so that it also preserves the energy. Since the time T depends only 
oil II/oIIlI; II/oIIlJi ||V'l|i°°7 ll/3||Li(-f ,f )j we can use the same arguments on the 
time intervals {T,2T), (2T, 3T), etc, to get a solution on (0, +00). Finally, in 
the case if 7^ 1, it is enough to rescale the time to get the right formulae, and 
to obtain the same conclusions. □ 



The following theorem claims the existence of a solution of (fTSj) in the non- 
cutoff case in some weak sense. 
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Theorem 2. Let /o S L\i^) be a nonnegative function. Let /3 satisfy the 
assumptions (HI) and (H2). Then, there exists a function g S L°°(M.^;C'^ D 
L°°(M)) which is a solution of I118\) in the weak sense defined by the forthcoming 
equation i2S^) . Moreover, g preserves the mass. 

Proof. We introduce, for n G N*, the cross-section 

MO) = I3{0) An = mmil3{e),n), 



and we denote by /" the solution of the problem ((18)) corresponding to the 
cross-section /3„ . This solution exists thanks to Theorem [T] For all n and all 
i > 0, fn{t) relies in L^, so that we can define its Fourier transform. Moreover, 
Qqb£;(/") also relies in L^. Hence, we can write the following equation: 

^I^^^ = r I /3„(^)/"/:(l + 5f~-{v') + SfHv'Me-"''^ ~ e'^^i)dvdvA9. 
at J^j^ Jr2 

We can split the integral on the right into three parts. The first part gives 
/3nWr/:(e-*"'« - e-^'^«)dz;d^;,d0 



pn{e) (r(ecos0)r(esin0) - r(or(o) 



dO. 



The second part can be evaluated using the inverse Fourier transform of the 
function /" which relies in L^ 



5_ 

2ti 



/3n(^)r/:rK)(e-'" « - e-'''^)dvdvAe 



MO) / f^f^Tif^){r^)e'v(vcose-v,sin( 



^^g-^^(^,cose-^,.si„e) _ ^~^v^\^ dr]dvdv^d9 



= ^ j\pn{o) / /"((e-^)cos0)r((e-?/)sin0)nr)('7)dr?d0 

I MO) I m -11 COS 0)f'\^sine)F{f^)(i^)df^d6. 

The third term can be computed in the same way. At the end we get that the 
Fourier transform of /" satisfies 
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dt 



J_[ Me) (r(ccos0)r(esin0) - r(or(o)) de 



S 



J'{n{v)dvde (22) 



_S_ 
27r 



/" COS0-T] sin 6*)/" (-^ sin 61 - r? cos 9) 



nrmdrtde. 



Note that here the mass is the quantity /(O). In addition 

Since the second moment of f^{v,t) is finite and conserved in time, its Fourier 
transform is two times differentiable, and satisfies 

sup|ia|^r(t,OllLoo < / v-'h{v)dv. 

t>o Jr 
Hence, we can use the Taylor formula at the order 2 

h{e) = h{o) + eh'{o) + [ (i - s)h"{se)ds 

Jo 

on the functions 9 ^ /" cos 6*) /" sin 61) , 61 /" ( - r?) cos 6*) /" ( - r?) sin 6*) , 
0^ ri^-ri cos 0)r{ri sin 9), 9^ cos9-r] sin 9)f"{-i sin 6* - cos 9) and 
^ f - V sin 9) f"{-r] cos 0). 

Using the notations m = J^fo{v)dv and e = J^v'^fo{v)dv, we get the fol- 
lowing estimates 



pn{0) (r(^cos^)r(^sin^) - r(or(o)) d^ 

< y ' 6i2^(6l)d6i (4|^pme + 2|^| 



(23) 



). 

y Pn{0) (r((e - ^) COS^)/" - 7?) Sine) 

- r(^ - r?cose)/"(77sine)V(r)(^)d^ 



< \i>{v)\m I 0^/3{0)d0 [4{\£, - t]\^ + \f]\^)me + 2(|^ - f]\ + \f]\)m^/^e^/^ 

(24) 
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and 



/" (C cose -J] sin 61)/" (-C sin 6 - t] cos 

- - 7? sin 0)r(-r; cos 0) )^(r)(7y)d0 



< \i^{v)\m / ' 02/3(0)d^ (4((|^| + \rj\f + \f]f)me + 2(|C| + 2\7^\)m^/^e^^^ 

(25) 

The right member of (j23p is integrable in time on any interval [ti, C M+, and 
the right-members of (j24l) and are integrable in (i,?]) on any [ti,t2] x 
with < ti < t2, since ip e C^(R). Therefore, to pass to the limit in it is 
enough for /" to converge pointwise on M+ x K. But inequalities ([Ml) and 
ensure that for all compact set /C C M, there exists a constant C depending 
only on /C, m, e, ■0! such that 

\r{ti,o-rit2,o\<c\h-t2\ vo<ii<i2, v^e/c. 

Then, thanks to Ascoli's theorem, there exists a function g E L°°{[ti,t2] x /C) 
such that, up to the extraction of a subsequence, 

II/" - 3llL'=([ti,t2]x/c) —-^ 0- 

All this being true for every ti, t2, JC, we deduce that g is well defined on M_|- x M 
and that g G L°°(M+; n C°(R)). We can therefore pass to the hmit in (1^ . 
Finally, we obtained the existence of a function g{^,t) which satisfies (|22p with 
the original cross-section /3, and such that 

g{t,0)^m Vi>0. 

□ 



4 Moments of the cutoff solutions, regularity of 
the non cutoff solution 

The second step of our analysis is the study of the regularity of the Fourier 
transform of the solution obtained in the previous Section. This can be done 
by investigating the moments of the solution to the cut-off equation. Since it is 
enough for our needs, we will limit ourselves to the fourth moment. Let 

.7r/2 

I5{e) cose sm^Bde 

J-n/2 

.7r/2 

A,= / /3(6')cos6i|sin6i|i+''d6i, - 1< 7? < 1, 

J-7r/2 
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J-7T/2 



'-7r/2 

and chose G N large enough such that, for ah n > np, 

l-Tv/2 2 

with /3n{9) = mm{(3{9), n). Our result is the following 

Theorem 3. Let j3{9) satisfy the assumptions (HI) and (H2), and let f"{t,v) 
be the solution of the problem il8\) with cross-section /?„ and with nonnegative 
initial datum fo G L\(R). Then there are explicit constants a>0, c>0,C>0 
such that we have 

sup lirWILj <(ll/o|Li+rf)e'^*, t>0, (26) 

n>n^ 

sup sup lirWILi < max{C, ||/o|Li}. (27) 

Moreover, a,c in i26]) depend only on ||/o||l1; ||/o||l1; ^"^^ whereas C 

in 1127}) depends only on ||/o||li> II/o||l1; '4' '^'^'^ ^r; md r/. 

Proof. Let us take (p{v) — as test function. Thanks to the symmetries of the 
kernel we obtain 

^ / v'f-{v,t)dv 

1 



2 , /3ni9) I {v'^ -~v^ + <^ - vt)rf:il + 6r\v') + 6f^iv:))dvdv,d9. 

From the collision rule (fTT|) it follows that 

(« ) + [v*) -V 

= («" + ^'*)(cos* 9 + sin^ 9-l) + 12v^vl sin^ 6'cos2 9 

+ Avv^ cos 9 sin 9{v'^ (sin^ 9 - cos^ 9) + vl (cos^ 9 - sin^ 9)) 

= -2 cos^ 9 sin^ 9[y'^ + v^) + 12i;2i;^ sin^ 9 cos^ 6i 
+ 4cos6'sin6'cos(26')ww*(i;J - v^). 

Consequently 

^iirwiLi 



7r/2 



-7r/2 
7r/2 



7r/2 
7r/2 



dvdvxd9 



(3r,{9) cos' 9 sin' l(vi + yj) f^^f- (l + 5r{v') + 5/"«)) 
+ 6 / ' (3 J9) cos' 9 sm'9 ff v'v'j'^ f: (l + S,r{v') + Sriv',)) dvdv,d9 

./-^/2 JjR2 ^ ^ 

/3„(6i) cos6'sin6icos(20) / / vvj"fl' (l + 5/"(w') + 5/"«)) d?;dv*d6' 



-7r/2 
/l +/2 +/3- 
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In what follows we always assume n > np. For the first and second terms 
/i, I2 we use 

r-iT/2 



1 f ' 

-A^ < / /3„ (9) cos^ 9 sin^ 9d0 < A* 

2 J-tt/2 



Ii<-\1 I l3{9)cos'^9sin^9d9 \ 2m I \v\^P{t,v)dv 

'-tt/2 I Jr (28) 



t/2 

and 

iirwiu- <m|iviiLo 

to get 

2 

= m2A,-mA||rWllLi, 

and 

/2 < 6(l + 2(5m||^||Loo)Ae2. (29) 
For the third term we compute by change of variable (v, v^,9) 1— > {v, w*, —9) 

13^6 /3„ (6*) cos 61 sin 6* cos(26') 

J-tt/2 

X - i;2)/"/^" (/"(^;') + /»«) _ _ dz;d«.d^ 

where 

v' = V cos 61 + sin 9, v'^ — —v sin 6' + cos 9. 

Recall 

v' = V cos 9 — sin 0, = f sin 6* + cos 6* 

and 

ii|;r(.,i)iu~<HiV''iiL~. 

It follows that 

!/"(«') + r«) - r(s') - ,r(s:)i < \~rw) - r(«')i + - rm 

<2m||V'||L-(|w| + |w,|)|sin0|. 

(30) 

The left hand side is also less than 2m||?/>||ioo (since / > 0). Thus 

ir(«')+r«)-r(^^')-r«)i < 2mm\L^f-'''{U'\\L^f-'''{w+\vA^)\ sm^r. 

(31) 

By and (|5T|) we obtain two estimates for /a: 

h<2m5W\\L^A ( ( \v\\v^,\{\v\ + \v^\){vl + v'')ri^Avdv^ 

J Jm? (32) 

and 

/3 < 2m{U\\L^Y-^{U'\\L^)^-nA^ 1 1 \v\\vA{\v\^ + \v.n{vl + z;2)/"/:Mi;d«* 

< K^.AmL^Y-'{u^\L^Y-'6A,{\\r{t)\\^^y-^ ■.= awrmLi)'^- 

(33) 
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Here if^.e > depends only on m and e. 

To prove the fist estimate ([2^ we omit the negative term Ii and use (P^ . 
([5^ and notice that A* < A to get 

^lirWLi < 6e'(l + 25mmL^)A + i^„,e<5||V/||L~ A|ir(t)||ii , < > 0. 
This gives by Gronwall lemma 

sup lirWLi < (||/o||Li+rf)e'^*, t>0 

where 

c = 6e2(l + 2Jm||V'||L-)A, a = K^JW\\l^A. 
To prove the second estimate (P71) we use (pS)) . and (1551) to see that 

^iirwiLi < C2 + cidirwL.)^ - mA4rm^., t > o. 

Here and below Ci, C2 and C3 depend only on m, e, V', ^, and 77. Ap- 

plying the following inequality 

Y"' <{l/e)^ +eY, a e (0, 1), F > 0, e > 



to 



gives 



and so 



Therefore, 



Y=\\fit)\\Ll, " = 



3 + J7 



^^\\rmLi<c3-^\\nt)\\Li. t>o. 



sup sup lirWILi < max <{ ^,11/0(011 



□ 



Passing to the limit n — >■ +00 in inequality (I27p we obtain 



Theorem 4. Lef (3 satisfy (HI) and (H2), and let g be the weak solution of 
hll^) defined in Theorem\^ with nonnegative initial data fo G L\. Then, g(t) is 
for all t and 

sup{||ff(i)lk~ + ||55gW|U^ + ||a|.5WI|L» + |l45WIU~ + l|a|45WllL~}<+oo. 
t>o 

Moreover, g conserves the energy, in the sense that dl^g{t, 0) = -e for all t > 0. 
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Proof. The conservation of the mass and inequality ^17} imply that there exists 
a constant C > which do not depend on n such that 



\dlJ"\\L^ 



1+ xl 



< c 

)<C 
)<C 

^ < c. 



Since L°°(M-|_ x M) is the dual space of the Banach space L^(M+ x R), the four 
sequences converge (up to the extraction of a subsequence) in L°°(R+ x R) 
weak-*; the limits can only be respectively d^g, d'^2g, d^^g and d'^^g (since the 



convergence in L'^ 



weak-* implies the convergence in the distributional 



sense). Moreover, we have the inequalities, for 1 < i < 4, 



\dhg\\ 



J+ xR) 



< 



liminf||5|.r|Uo 



Finally, we have the embedding 



so that g e L°°{R+;C^{R)). 

It remains to prove that the energy is conserved. Let us fix some time to > 0. 
It is clear that for all integer n, we have 



\\dir{to,.)h 



< e. 



Therefore, up to the extraction of a subsequence, there exists a function h S 
L°°(R) such that 



weak — * L° 



But it is clear that h = d'^2g{to,.), since for all function (p ^ C° 
compact support, 

: Jr 
or, after two integrations per part on both sides, 

4r(io,e)0(Ode^ / 4.9(io,e)0(Ode 



with 



Let us define an approximation of the Dirac measure 
We have 



P 




otherwise 



4r(*o,0)-45(to,0) 



< 



4r(io,o)- / <i>p(c)4r(to,c)d^ 



pio (a|2r(io,0"a|.<?(io,o)de 



<5'p(O45(io,Od$-9|25(to,0) 
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The first and the third terms converge toward when p converges to infinity, 
independently of n since 

ii4/"(*o,.)iu- <c 

with C > independent of n. As for the second term, once p has been fixed, it 
converges to since <i>p e L^(R), and the result follows. □ 

5 The grazing collision limit 

We are now in a position to perform the grazing collision limit in equation 
(|18p . We can now make precise assumptions on the asymptotics of the grazing 
collisions, namely in letting the kernel /3 concentrate on the singularity = Q. 
We will introduce a family of kernels {/3e(|6'|)}£>o satisfying hypotheses (HI) 
and (H2), with 

lim / pMe\)e'^de = i (m) 

e^o+ J a 

This can be obtained in several ways, for example taking, for < /i < 1 

^^(1^1) = 0<|^|<eva-.), 

^-(I^D^^j^p^ elsewhere. 

Let 5£ be the weak solution of the problem (fT8)) in the sense that it satisfies 
equation f^^ . where P{9) has been replaced by (3^{9). 

Theorem 5. Let ji{9) satisfy assumptions (HI), (H2), and let /3e(0) satisfy 
\34-^ . Let be the weak solution of the problem il8\} where /3(0) has been 
replaced by l3^{0), with the nonnegative initial data fo satisfying /p G L4(IR). 

Then, for all T > 0, there exists a function g G L°°{0,T;W^'°°{R)) such 
that, up to the extraction of a subsequence, 

- 5lU-([ii,t2]x/c) ^ yO<ti<t2<T andlCcM. compact, 



dlge dLg in L°°((0,T)xR) weak-*, 1 < i < 4, 
and such that g is a solution of the Fourier form of the equation 



[ h{v)dv) ^{vh{l + 5h)) +( [ v'h{v){l + 6hiv))dv) P^, (35) 



dt J dv \Jg^ J dv 

Remark 6. If (5 = 0, equation (|55|) reduces to the classical linear Fokker-Planck 
equation. 

Proof. To pass to the limit, we need some regularity on the solution t^g. However, 
Theorem S] is not sufficient, since the bound on the fourth derivative of 
depends on wich becomes an unbounded quantity when replacing /3 by /?£ 
and as £ — 0. Therefore, we use inequality (j26p to obtain a regularity result on 
g^ which is the same than Theorem 01 except that it works only on any time 
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interval (0,T). The gain is that the constants used remain bounded as e — > 
when replacing /? by 

Let us fix some T > 0; 5^ is four times differentiable for almost any time in 
(0,T), and these derivatives are bounded uniformly in time and independently 
of e, provided the initial datum /g relies in L\{M.). We now act as if (3{6) was 
integrable, and then we will obtain the result by an approximation argument. 
A Taylor expansion in e under the integral sign gives 



dt 
1 



m 



+ {{V - 0' - rf)9e{i - v)9M9e{0) + Ui09eiv)9ei-V) 

+ e9':i-v)9eM9eiO+^vUi09eiv)9'ei-v))dri^ +9^0{e^) dO. 

(36) 

Since (3(9) is an even function, the first-order terms vanish. 

By Theorem |3] there exists a constant At > which do not depend on e such 
that 

sup \j2\\dh9s{t,.)\\L^ \ <Xt. 

Using equation p6p . we see that the family (ge)^ is equicontinuous, so that we 
can use Ascoli's theorem, which says that there exists a function g G L°°{{0, T) x 
R) such that, up to the extraction of a subsequence, 

for all < ti < ^2 < 7" and all compact set /C C M. In addition, all the 
results of Theorem H] are still valid for g on (0,r). Therefore, thanks to both 
the uniform convergence for g and the convergence in L°°((0,T) x R) weak-* 
for its derivatives, we can pass to the limit in equation (1361) . and we get, using 
classical formulae on the Fourier transform and the conservation of mass and 
energy, that g satisfies the equation which is the Fourier transform of equation 

(133). □ 



6 Other results 

To give a sense to the collision kernel, we have been forced to mollify a part of 
our original equation 

^ = QQBEi.f){t,v), teR+, veR, ^^^^ 
f{0,v)^fo{v), 



18 



where 



QQBE{m,v) ^ m {ff'M+sm + su) 

However, assuming that a solution to this equation exists and is regular enough, 
we are able to obtain some interesting consequences. 

6.1 if-theorem 

The first one is the iJ-theorem. Under the right assumptions on the solution / 
of equation p7p , we can define 

H{f)^ £(^^{1 + Sf) log(l + 5f)-f log d^; 

and 



Dif) = \ l_[ m j^^ r(//:(i + 5f){i+5u), 

fU{l + 6f){l + 6fS)'^vdvAe 



with 



(a - 6) log(a/6), a > 0, 6 > 0; 
T{a,h) = {+00, a>0, 6 = or a = 0, 6 > 0; (39) 

0, a = 6 = 0. 

Then, / satisfies the so-called 7J-theorem: 

Theorem 7. Let f he a solution of the problem i37\} . with fo £ L\ogL{M.), and 
assume that H{f) and D(f) are well defined. Then 

H{f{t, .)) = H{fo) + f D{f{s, .))ds yt > 0. 
Jo 

Consequently, the entropy H is increasing along the solution. 

6.2 Regularity of the solution 

Using the i7-theorem, we can give an a priori estimate on the solution / of p7p . 

Theorem 8. Let /3 satisfy the properties (HI), (H2) and (H3), and let f{t,v) be 
a solution of the problem 1137]) with the initial datum fo G L\ flLlogL. Assume 
that f satisfies the H -theorem. Then we have 

iog(i + ^/)eLL(K+;^''/' 

If in addition f G L°°(M+ x K), then we have 

fGLlMv.H'^'^ 
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Proof. We will do the computations as if the cross-section and the function / 
were smooth. Using the classical changes of variable {v,v^,,d) i— ?> {v',vl,—6) and 
{v,v^) I— >■ which have unit jacobian, we have: 



1 /"^ 

2 



^ , f'mi + 5f){i+sf,) ^ 



/5W / (/'/:(! + -J/)!! + 'J/*) - fMl + Sf'){l + Sfi)) 
X log (//.(I + <5/')(l + <5/:))dz;di;,d0 

m I (/'/:(! + + 5f*) - fMl + smi + sf:)) 

X log (/(l + (5/'))dwdv,d6' 

- /(I + 5f) + /'(I + <5/))dwdi;*d0 

+ j\ m j^^ + -5/:) (/(I + -^/o ~ /'(I + 5f)yvdvAe 

= /l+/2. 

The term I2 can be treated easily, because it can be written as 

m I Mi+sf:Kf~f )dvdv,d9, 



and the presence of / — /' involves strong cancellations. In fact, the term I2 
verifies 

/.HI/llL/%w(i-^)d.. (40) 

To prove (PO]) . consider that 

m I f.u - ndvdvAe+6 r m ( m'm - ndvdv^de 



Thanks to the change of variable {v,v^,,9) 1— {v' ,vl,—9), we see that /| = 0. 
On the second part of I2 , we use the change of variable v t-^ v' with and 9 
fixed, which jacobian is 

dv _ 1 
dv' cos 9 ' 
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Therefore 

and from this (j40D follows. 
Using now the inequality 

a; log x + y> {y/x - y/y) , Vx, y > 0, 

y 



we obtain 



h> [[m I + -5/:) + V) - + <5/))'dz;clz;,cl0. 



Therefore 



/3(0) / /,(l+<5/:)(V/(l + <5/')- + 5f)) AvdvAe < D{^)+Cl\\f\\l^ 

(41) 



where 

ci = /' /3(6') f ^— - 1 ) d6l > 0. 



008 6* 



Now, we write 



For < a < X. let 



, , _ \/a;(l + (5a) — \/ a(l + ^x) 
~ log(l + <5a) - log(l + bx) ■ 

Then, there exists some constant C2 > that does not depend on x or a such 
that 

\<^{x)\ > C2 > a. 

>From this inequality we deduce that 

C2 /'^m / /4l+'5/:)(l0g(l + 5/')-l0g(l+5/))'<i?(/)+Cl||/||il 

It has been shown in \V\ that if F is a real function such that F(/) e L^(M) 
satisfying 

m\^^ f*{Fif') ~ Fif))' < Dif) + cWfWl, 
then the following inequality holds: 

mf)\\%.,.<D{f) + c^\\f\\i.. 

Taking F{f) = C2 log(l + Sf) the result follows □ 
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6.3 Moments of the solution 



To study the grazing collision limit, we needed some regularity on the Fourier 
transform of the solution. This is equivalent to have a uniform bound on some 
higher moment of the solution. In proving Theorem |21 we used the regularity 
of the mollified part, more precisely the fact that this part was in C^. In fact it 
could be enough to use the regularity which follows from the i?-theorem. 
Indeed, the terms that raise problems in the proof of Theorem [3] are Ii and l2- 
Let us see how to treat the first one. We have 



I3{e) sin6'cos^ 



vv^v^ - vDfMfiv') - f{v'j)dvdv,de 



/' 



/3(6') sin 6* cos^ 61 / vv^{v'^ ~ vl)f{v)f{v^) 



v{l ~ cost 



d?;du,d6' 



[v{l - cos 6*)]" 
< / /3(6')|sin6i||l -cos6i|"|cos3 



-dvdVi 



i/« 



d0. 



In order to deal with moments not exceeding the fourth one, we use Holder's 
inequality, with p such that (3 + a)p = 4. Consequently q = Moreover, 
in order to recognize the semi-norm of / in the Sobolev space H'^^^ in the last 
term of the product, we need to set 



and thus 



aq — 1 + V, 

l + v 
a = 



(note that < a < 1). With these constants, we have 



\fW)~f{i')\' 
\v{l - cos^)!"? 



dvdvi, 



1/9 



1 



< 



11-008 6*11/9 

1 2 



i/p 



'/(f)-^(:?^ .d(i-cos.K 

\v{l - cos6l)|"'? ^ ' 



|1 — cos 6*1 1 1 sin 6*1 « 



rll/WILoo'l/WI^. 



/2- 



Finally 

\h\<c [ m\cos^' 



I sin 6*1^ 



|l-cos6ir"id6l 



v^f{t,v)dv ||/WIU~I/WI^. 



The integral in 6 is finite if and only if 

2^-1 
V 9 

Recalling that 



- - - > -1. 

q 



q = 



1 - a 
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and that 

l + y 

a = , 

5 + 1/ 

this request is equivalent to 

+ - 4 < 0, 

which is verified for —1 — a/S < v < \/5 — 1. Combining this with the previous 
constraints on tire parameter i/, we obtain that our bound of the fourth moment 
works for all v verifying 

Kv <Vl-l. 

The treatment of I2 is quite more simple, since it requires only some changes 
of variable. Indeed, using the changes of variable w 1— > w — sm9 and 9 1— > —9 



13 {9) sin 9 cos'^ 9 / vv4v^-vi)f{v)f{v^){f{d')-f{v))dvdvA9 
;3(6l)sin6lcos3 6' / vv4v^-vl)f{v)f{v^)f{i')dvdv^d9 
;3(6l) sin 61 cos^ei- 

{v + sin6')t;,((v + sin9)^ — v1)f{v + sn\9) f {v^) f {v)dvdv^:d9 

I3{9) sin6'cos^ 9- 
{v - sm9)v^{{v ~ sin6')2 - vl)f{v')f{v^)f{v)dvdv^d9 
This implies 

I2 < C(m,e) (^J' I3{9)sin^0d9^ + to^ (^J v^f{t,v)dv 

We proved 

Theorem 9. Assume that the cross-section /3 satisfies (HI), (H2), with 1 < 
v < 1. Assume that there exists a solution f to the problem iS7\) . Assume 

that f is regular, in the sense that \\f{t)\\L°°(M) and \\f{t)\\^^ are in _LP(R+) for 
some p big enough; then there exists some constant A > such that f satisfies 



<max{A,||/oLi} 



7 Conclusions 

In this paper we investigated the asymptotic equivalence between the (mollified) 
Kac caricature od a Bose-Einstein gas and a nonlinear Fokker-Planck type equa- 
tion in the so-called grazing collision limit. The limit equation differs from the 
analogous one present in the literature [H] , since in our case the linear diffusion 
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has a diffusivity which depends on the solution itself, in order to guarantee the 
conservation of energy. Our analysis refers to a mollified version of the equa- 
tion, due to the difficulties of handle the third order nonlinearity present in the 
Bose-Einstein correction. A further inside on the true model, done in the last 
part of the paper, shows that a proof of the boundedness of the solution would 
be sufficient to avoid the presence of the moUifier. 
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